Tailoring Magnetism in Quantum Dots 
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We study magnetism in magnetically doped quantum dots as a function of confining potential, 
particle numbers, temperature, and strength of Coulomb interactions. We explore possibility of 
tailoring magnetism by controlling the electron-electron Coulomb interaction, without changing the 
number of particles. The interplay of strong Coulomb interactions and quantum confinement leads 
to enhanced inhomogeneous magnetization which persist at higher temperatures than in the non- 
interacting case. The temperature of the onset of magnetization can be controlled by changing the 
number of particles as well as by modifying the quantum confinement and the strength of Coulomb 
interactions. We predict a series of electronic spin transitions which arise from the competition 
between the many-body gap and magnetic thermal fluctuations. 

PACS numbers: 75.75.+a,75.50.Pp,85.75.-d 



Magnetic doping of semiconductor quantum dots 
(QDs) provides an interesting interplay of interaction ef- 
fects in confined geometries [3, 0, M 13, H, IE 0, Hi and 
potential spintronic applications [9j. In the bulk-like di- 
lute magnetic semiconductors the carrier-mediated fer- 
romagnetism can be photoinduced IIOL lll | and electri- 
cally controlled by gate electrodes QJI], suggesting pos- 
sible nonvolatile devices with tunable optical, electrical, 
and magnetic properties QDs allow for a versatile 
control of the number of carriers, spin, and the effects 
of quantum confinement which could lead to improved 
optical, transport, and magnetic properties as compared 
to their bulk counterparts [l|, |l3|, [l4[ ■ Unlike in the bulk 



structures, adding a single carrier in a magnetic QD can 
have important ramifications. An extra carrier can both 
strongly change the total carrier spin and the temper- 
ature of the onset of magnetization which we show can 
be further controlled by modifying the quantum confine- 
ment and the strength of Coulomb interactions. 

We study the magnetic ordering of carrier spin and 
magnetic impurities in (II,Mn)VI QDs identified as a ver- 
satile system to demonstrate interplay of quan tum con- 
finement and magnetism 

US BE [3 EMI- Because 



Mn is isoelectronic with group-II elements it does not 
change the number of carriers which in QDs are con- 
trolled by either chemical doping or by external electro- 
static potential applied to the metallic gates. The latter 
allows confinement of the carriers in a dot with tunable 
size and shape • By using real space finite- temperature 
local spin density approximation (LSDA) fl9| we study 
temperature (T) evolution of magnetic properties of QDs 
over a large parameter space. This approach allow us to 
consider QDs with varying number of interacting elec- 
trons (N) and Mn impurities (N m ) which already for 
small N and N m becomes computationall y in accessible to 
the exact diagonalization techniques HI, 2(| • We extend 
the previous studies of Coulomb interactions in magnetic 
QDs with N m = 1, 2 at T = [H| and T > results us- 



ing either Thomas-Fermi approximation or by applying 
Hund's rule with up to 6 carriers [13] • We reveal that the 
interplay of strong Coulomb interactions and quantum 
confinement leads to enhanced inhomogeneous magneti- 
zation which persist at higher temperatures than in the 
non-interacting case and the bulk structures ID, 1Z|- We 
refer to such a spin-polarized state in QD at zer o ap plied 



magnetic field as "ferromagnetic" state [la, Il3, ll 



Here we focus on magnetic QD in zero applied 
magnetic field described by the Hamiltonian H = 



He 



H m + H ex , with the electron contribution H e 



—r , where U is the 

V 3 

is the electron effective mass, and 



Planck constant, m 
Uqd(?) is the confining potential of a three-dimensional 
QD. The last term in the equation is the repulsive 
electron-electron (e-e) Coulomb interaction screened by 
the dielectric constant e, — e is electron charge, and 7 ac- 
counts for reduction of Coulomb strength due to screen- 
ing effects of the gate electrodes 21| . The Mn Hamilto- 
nian is H m =Yli v JfpMj-Mjt, where J AF is the direct 
Mn-Mn antiferromagnetic coupling. The z-component of 
Mj of impurity spin satisfies M z = — M, — M + 1, . . . , M, 
where we choose z as the quantization axis and M = 5/2 
for Mn. The electron-Mn (e-Mn) exchange Hamilto- 
nian is H ex = -J s dJ2ijSi ' M iH r i - where J sd 
is the exchange coupling between electron spin Si, at 
Yi = (pi,Zi), and impurity spin Mj, at Rj. An effec- 
tive mean field Hamiltonian describing electrons can be 
obtained by replacing the Mn spins, that are randomly 
distributed, with a classical continuous field H® s = 
H e — J2 i Jsdn m ^(M z (rj)), where n m is the averaged 
density of Mn, and a = ±1 for spin up (|), and down 
(|). The effective magnetic field seen by the electrons 
is the mean field induced by Mn. Assuming that im- 
purities are in equilibrium with thermal bath it follows 
(M z {vi)) = MB M (Mb(Ti)/k B T) where B M (x) is the 
Brillouin function [22j |. ks is the Boltzmann constant. 
Here b{v t ) = -Z Mn J AF (M,(r<)) + J^K^) - n x (r 4 )]/2 



2 



is the effective fieid seen by the Mn 23[. The first term 
in 6(r,*) describes the mean fieid of the direct Mn-Mn an- 
tiferromagnetic coupling [l6j]. Z Mn is the averaged Mn 
coordination number, and n (T {T i ) is spin- resolved elec- 
tron density. We decompose the planar and perpendicu- 
lar components of the confining potential of a single QD, 
and fit it to a realistic QD potential 24|. The result- 
ing potential, Uqd, is a sum of a two-dimensional (2D) 
Gaussian Vqd = Vo exp(— p 2 /A 2 ) and one-dimensional 
parabolic potential Vq D = m*Vl 2 z 2 /2, where p = (x,y). 
For Vqd we find that Gaussian potential is more realistic 
than usually studied parabolic potential. Here Vq and 
are the planar depth of the QD minimum, and the char- 
acteristic subband energy associated with the perpendic- 
ular confinement. In typical disk-shaped QDs, and low 
density of electrons, only the first subband is filled. Af- 
ter expanding the QD wave functions in terms of its pla- 
nar ipi a (p) and subband wave function we project 
-ff| ff into a two-dimensional Hamiltonian by integrating 
out £(z). In LSDA the two-body Coulomb interaction 
can be written as sum of Hartree potential Vh and spin 
dependent exchange-correlation potential V£ c . We use 
Vosko-Wilk-Nusair exchange-correlation functional [lj|, 
and express the Kohn-Sham Hamiltonian as 

Hks = + Vqd + jV H + iV° xc - ^h sd (p), (1) 

where 



h s d{p) = Jem / dz\H(z)\ B M 



,a D fMb(p,z) 



(2) 



and J em — J sc in m M is the e-Mn exchange coupling. 
The Kohn-Sham eigenvectors and eigenvalues of Eq. |T]), 
ipnaip), and e ncr are calculated numerically. 

We illustrate our findings and the iterative solutions 
of Eq. (P) for (Cd,Mn)Te QD. The material parameters 
are J sd = 0.015 eV nm 3 , m* = 0.106, e = 10.6 [3, and 
we choose Zmu J AF = 0.02 meV. The planar (x,y), and 
perpendicular (z), dimensions of the QD are taken as 42 
nm and 1 nm with n m = 0, 0.025, 0.1 ran -3 . In the 
central region of QD of area ia*g, where a* B — 5.29 nm 
is the effective Bohr radius in CdTe, n m = 0.1 nm~ 3 cor- 
responds to w 10 Mn atoms. For a planar confinement, 
Vqd, we consider a Gaussian potential with Vq — —128 
meV, and A = 38.4 meV, corresponding to ujq — 27 meV. 
Here u>o is calculated by expanding Vqd in the vicinity of 
the minimum which yields Vqd = Vo + m*uj 2 p 2 /2 + . . ., 
with the strength o> = ^/2|Vo|/m*/A. 

In QDs electron density is inhomogeneous, implying 
that both the electron spin density, n^{p) — n^{p), and 
Mn-magnetization density (M z (p)) = Mh sc [(p)/J em are 



inhomogeneous. For N = 8 and 



0.1 nm , we 



show the self-consistent spin density in Fig. [TJa), and 
Mn-magnetization density in Fig.QJb). Outside the QD, 
nicrip) decays exponentially, and an effective field b(p) 




M 




FIG. 1: The spatial density profile of electron spin density 
(a) and Mn-magnetization (b) for QD with N = 8 at T = IK. 
Coordinates (a;, y) axe expressed in effective Bohr radius. 



seen by Mn becomes negligible. This is consistent with 
vanishing (M z (p)} at the QD boundary [Fig. QJb)]. 

We next turn to spatially-averaged quantities such 
as Mn-magnetization per unit area A, (M z ) = 

J d 2 p(M z (p)}, electron (spin) polarization P = (iVf — 
N\)/N, and the z-component of the total spin of elec- 
trons, s z = (N T - iVj/2. In Fig. Owe show (M z ) as 
a function of TV for n m = 0.025 nm~ 3 , and both non- 
interacting (7 = 0), and interacting (7 = 1) electrons. 
The magnetic behavior of QD can be described based on 
the interplay of the many-body spectrum (determined 
by the shell structure of the confining potential, and e-e 
Coulomb interaction) and the strength of e-Mn exchange 
coupling J em . In the following we summarize the spin 
structure of the QD in the absence and presence of J em 
with 7 = 0, 1. 

i) Jem = 0: The shell structure of the 2D Gaussian 
potential is shown in Fig.[2ja). The energy gap between 
s-, p-, and d-orbitals is characterized by u>q. In contrast 
to 2D parabolic potential [H, Hj], ci-shell levels are not 
completely degenerate, and therefore we focus on N > 6 
states. Degenerate levels d + and d- are separated by an 
energy gap (1.5 meV) from do -level, where ±,0 refer to 
angular momentum l z = ±1,0. However, e-e interaction 
changes the structure of (i-shell as it overturns the order- 
ing of the d-orbitals, e.g., the Kohn-Sham energies of pair 
of degenerate d+, and cL are below do (with energy gap 
w 1 meV). Because of <i-shell overturning, N = 10,12 
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FIG. 2: The averaged magnetization per unit area (M z ) as a 
function of number of electrons NatT= IK and Mn-density 
n m = 0.025 nm -3 for non-interacting (7 = 0, empty trian- 
gles) and interacting (7 = 1, filled triangles) electrons. The 
ground state of the QD switches between ferromagnetic and 
antiferromagnetic states as function of N. Coulomb interac- 
tion changes the state of N = 8 from antiferromagnetic to 
ferromagnetic state, (a) The schematic single particle levels 
of a 2D Gaussian confining potential, (b) The z-component 
of the total spin of electrons, s z , as a function of N for 7=1. 
(c) Antiferromagnetic-ferromagnetic transitions for N = 8 as 
function of n m and 7. 
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FIG. 3: Temperature evolution of Mn-magnetization per unit 
area (M z ) (a), the electron polarization P (b), and the free 
energy difference AF between ferromagnetic and antiferro- 
magnetic QD (c). At low T, N = 4, 8 form half-filled shells 
with P = 2/4, 2/8. T — T* , characterizes vanishing of (M z ), 
P, and AF. 



form closed shells with s z — 0, and N — 7,9,11 form 
open shells with s z = 1/2. The N — 8 corresponds to 
half- filled shell with s z — 1, and electron polarization, 
P = 2/8. The evolution of s z as a function of N for 
7 = 1 is shown in Fig. (2fb) (empty circles). 

ii) Jem 7^ 0: 7 = and increasing e-Mn coupling to 
Jem = 3.75 meV, leads to transitions P = 0/8 — * 2/8, 
and P = 1/9 -> 3/9 for N = 8, and 9, whereas N = 10 
shows P = 2/10. For 7 = 1, the dependence of P on J em 
is negligible in s-, and p- shells. In contrast, in d-shell, we 
find transitions P = 1/9 — * 3/9, and P = 0/10 -> 2/10 
for N — 9 and N = 10 at low T. However, we find 
no change in P for N = 7,8,11, and 12. Figure [2Jb) 
(filled circles) show s z as a function of for n m = 0.025 
nm~ 3 (J em = 0.94 meV). Increasing the density of Mn 
to n m — 0.1 nm~ 3 (J em — 3.75 meV) does not change 
s z . 

In Fig. [21 we observe that (M z ) = in closed shells 
(iV = 2, 6, 12) for both 7 = and 7=1, because of well 
separated s-,p-, and d- orbitals due to large cuo(~ 30J em ). 
Comparing (M z ) between 7 = and 7=1 one can ob- 
serve that the e-e interaction stabilizes the ferromagnetic 
state due to the spin Hund's rule. This condition is easily 
satisfied for open shells where the maximum electron po- 
larization is obtained in half-filled shell with A^ = 4. The 



N = 8 state (recall Fig. [T]) is more interesting. At 7 = 
electrons fill single particle levels following the Pauli ex- 
clusion principle. Even with J em — 0.94 meV (smaller 
than single particle e-h excitation gap), N = 8 forms 
closed shell and P = (M z ) = 0. In the case of 7 = 1, 
and because of (i-shell overturning, polarized electrons 
in d + and gL give P = 2/8, and finite (M z ). We also 
see that the maximum (M z ) occurs at A^ = 9 because 
Jem = 0.94 meV induces three polarized electrons in d- 
levels. Figure [2tc) reveals the dependence of magnetic 
transitions on 7, and n m . With increasing ZmuJ AF , the 
transition to ferromagnetic state occurs at larger 7. Our 
findings clearly demonstrate that the magnetism induced 
by strong Coulomb interaction can be controlled by the 
electric gates or by changing the semiconductor host (and 
thus changing e) without changing the number of carri- 
ers confined in QD. We also suggest that because of the 
sensitivity of the (M z ) to the electronic spin transitions, 
the former can be used to infer the spin of electrons, 
and could be potentially applied to manipulation of spin 
qubits in semiconductor nanostructures [9j . 

We next examine the temperature dependence of mag- 
netism in QDs. In Fig. [3l we show (M z ) (a), P (b), and 
the free energy difference AF between ferromagnetic and 
antiferromagnetic states (c) for N = 4 and N = 8. The 
suppression of (M z ), shown in Fig. [3Ja) is accompanied 
by a series of spin transitions in electronic states and sup- 
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FIG. 4: T* as a function of N for interacting electrons 7 = 1, 
ljo = 27 meV, and Vo = —125 meV. Inset: The dependence of 
T* on coo for N = 1, N = 4 with 7 = 0, and 7=1. w = 27 
meV is marked as a dotted line. There is an optimal confining 
potential which maximizes T* . 



peratures where no such ordering exist in the host ma- 
terial [26]. In the limit of small u>o, where the Coulomb 
interaction among particles is the largest characteristic 
parameter of the QDs, we find magnetism substantially 
stronger than predicted from the non-interacting picture. 
In contrast to the carrier-controlled ferromagnetism in 
the bulk-like structures l3, 11, 12 1, we reveal that mag- 
netism in QDs, can be tuned even at the fixed number of 
carriers by gate voltage which controls the inter-particle 
Coulomb interaction screening. For potential spintronic 
applications based on II- VI magnetic QDs, we antici- 
pate that it is possible to further increase the magne- 
tization and the temperature at which it vanishes. In 
addition to exploring a larger hole-Mn exchange cou- 
pling in (II,Mn)VI QDs, it would also be advanta- 
geous to consider (II,Cr)VI QDs as there is a support 
for the room-temperature ferromagnetism in their bulk 
counterparts 
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pression of P. At low T, the spin triplet is realized as 
the ground state of the N = 4, and N — 8 open p-, and 
d-shells (due to Hund's rule). We define a characteristic 
temperature, T* , at which (M z ) = P = AF — 0. 

In Fig. H we plot T*(N) for uj = 27 meV and 7 = 1 
which decreases non-monotonically with N. The inset 
shows T*(lo q ) for N = 1, and N = 4 (7 = 0, 1). At low 
u)q, the e-e interaction strongly enhances T* , while at 
large loq the effect of confinement potential is dominant. 
Thus we find T*( 7 = 1) -> T*( 7 = 0) with increas- 
ing ujq, which in turn gives rise to a peak in T* . Sev- 
eral trends in calculated T*(N, ojq) can be obtained from 
a perturbative approach by approximating 2D Gaussian 
with 2D parabolic potential. Near (M z ) = P = for 
QD with one valence electron in s-, p-, or d-shells, we 

find T* = J emv /g^[/d 3 r|^(r)| 4 ] 1/2 , where ^ is 
the wave-function of the highest occupied orbital, and 
J AF = 0. For a given uio, T* decreases with N, e.g., 
T^ =3 = 0.7T^ =1 and T^ =7 = 0.6T^ =1 . One can also 
show that T* oc y/ujo, consistent with bound magnetic 
polarons [251 ]. 

In conclusion we have investigated the existence of 
magnetism in magnetically doped QDs, as function 
of particle numbers, confining potential, temperature 
and strength of Coulomb interactions, using finite- 
temperature LSDA. Our results show that QDs embed- 
ded in magnetic semiconductor host can be considered 
as ferromagnetic centers which exhibit spatial ordering 
in spin density and magnetization, even at elevated tem- 
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